Muon Spin Rotation Measurement of the Magnetic Field Penetration Depth in 
Ba(Fe .93Co 07)2As 2 : Evidence for Multiple Superconducting Gaps 
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We have performed transverse field muon spin rotation measurements of single crystals of 
Ba(Feo.93Coo.o7)2As2 with the applied magnetic field along the c direction. Fourier transforms 
of the measured spectra reveal an anisotropic lineshape characteristic of an Abrikosov vortex lat- 
tice. We have fit the /xSR spectra to a microscopic model in terms of the penetration depth A and 
the Ginzburg-Landau parameter k. We find that as a function of temperature, the penetration 
depth varies more rapidly than in standard weak coupled BCS theory. For this reason we first fit 
the temperature dependence to a power law where the power varies from 1.6 to 2.2 as the field 
changes from 200G to 1000G. Due to the surprisingly strong field dependence of the power and the 
superfluid density we proceeded to fit the temperature dependence to a two gap model, where the 
size of the two gaps is field independent. From this model, we obtained gaps of 2Aj = 3.7A;_bT c and 
2A2 = 1.6A:_bT c , corresponding to roughly 6 meV and 3 meV respectively. 

PACS numbers: 73.43. Nq, 75.30. Kz, 75.30.Sg, 75.40.Cx 



Recently, a new family of high temperature supercon- 
ductors were discovered p] based on layers of FeAs. Sev- 
eral groups of these structures have been found, includ- 
ing the so-called 1111 compounds, such as LaFeAsO, and 
the 122 compounds, such as BaFe2As2- Superconductiv- 
ity was discovered in these systems by hole-doping, such 
as partially replacing oxygen by fluorine in LaFeAsO or 
partially substituting potassium for barium in BaFe2As2- 
The resulting superconductors have remarkably high 
T c 's, for example, up to 56K for LaFeAsfpi-zF^HJSIl] 
and up to 38K for (Bai_ x K x )Fe 2 As2|5J EJ 0. 

At room temperature the parent compounds have 
a tetragonal structure, which then undergoes an or- 
thorhombic distortion somewhat below 200 K. Either 
at or just below the structural phase transition, there 
is a magnetic phase transition to an antiferromagnetic 
(AF) ground state. With doping, the AF transition is 
supressed, and superconductivity emerges. The proxim- 
ity to AF magnetic order has led to a belief that Fe 
spin fluctuations are important for developing the su- 
perconducting ground state |5]. More recently, it was 
found that superconductivity can be induced in these 
systems through electron-doping through substitution of 
iron in the FeAs layers; both La(Fei_ a; Co a: )AsO[9] and 
Ba(Fei_ x Co a; )2As2[Tni[II] become superconducting with 
Tc's as high as 14 and 23 K respectively. In contrast to 
the cuprates, the iron pnictides are apparently surpris- 
ingly robust against in-plane disorder. 

At this point there is no consensus on the pairing sym- 
metry in pnictides. In the 1111 system most experiments 
indicate a fully gapped Fermi surface, although the pres- 



ence of magnetic rare earth atoms and the lack of large 
single crystals complicates the interpretation [12]. In the 
122 systems the situation is also unclear. Tunnel diode 
resonator measurements [T3J[H] of Ba(Fei_ :c Co :r )2As2 ex- 
hibit power law temperature dependences for the pen- 
etration depth AA(T) ~ T" with n between 2 and 
2.5 which have been interpreted in terms of gap nodes. 
Similar measurements of Bai_ a; K 2 ;Fe2As2[15 exhibited a 
power law n of about 2 for low (T/T c < 0.25 and were fit 
with a two gap s-wave model with somewhat unphysical 
parameters. Muon spin rotation measurements of sin- 
gle crystal Bai_ a: K a ;Fe2As2 exhibited phase separation 
into magnetic and superconducting regions |16| 117] , with 
different temperature dependences for the superfluid den- 
sity (n s oc 1/A 2 ) depending on the value of x. NMR 1/Ti 
measurements [TB] show no Hebel-Slichter peak below T c 
which generally indicates non-s-wave pairing, whereas 
the 59 Co and 75 As Knight shifts decrease below T c for 
fields both along the c-axis and in the a-b plane; behavior 
which is consistent with s-wave pairing. 

Muon spin rotation (//SR) is a powerful local micro- 
scopic tool for characterizing the magnetic properties of 
materials, in superconducting or other states. A thor- 
ough description of the application of /iSR to studies 
of superconductivity can be found elsewhere[19 . In a 
transverse field (TF) /iSR experiment, spin polarized pos- 
itive muons are implanted one at a time into a sample. 
The muon spins precess around the local magnetic field, 
and decay into a positron, which is preferentially ejected 
along the direction of the muon spin at the time of de- 
cay (as well as two neutrinos which are not detected). 
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In the presence of a vortex lattice, the spatial variation 
of the magnetic field distribution results in a dephasing 
of the muon spin polarization and a relaxation of the 
precession signal. A Fourier transform of the spin polar- 
ization function essentially reveals the field distribution 
which exhibits a characteristic Abrikosov lincshapc. The 
lincshape (or equivalcntly the relaxation function in the 
time domain) depends on the lattice geometry, magnetic 
field penetration depth A, coherence length, £, and the 
amount of lattice disorder. As a result, careful analysis of 
the relaxation function allows these microscopic parame- 
ters to be determined in the vortex state. Such measure- 
ments demonstrated the presence of gap nodes character- 
istic of d-wave pairing in high quality single crystals of 
YBa 2 Cu30 6 .97[Tn]. In ceramic samples this anisotropic 
lineshape is generally not observed, rather the broadened 
line is generally well described by a gaussian distribution; 
however, the width of this distribution (the Gaussian re- 
laxation rate) a has been shown to be proportional to the 
superfluid density a cx n s cx l /A 2 |20ll2Tj . Previous stud- 
ies of cuprates found that extrinsic effects in ceramics 
can result in the correct temperature dependence of the 
superfluid density being masked; for this reason, reliable 
measurements of the superfluid density require the use 
of single crystals and the observation of an anisotropic 
lincshape characteristic of a vortex lattice. 

A single crystal of Ba(Feo.93Coo.o7)2As2 was grown 
from self flux as described in detail elsewhere [11 . This 
crystal of roughly 1cm 2 area was mounted in a helium 
gas flow cryostat on the M20 surface muon beamlinc 
at TRIUMF, using a low background arrangement such 
that only positrons originating from muons landing in the 
specimens were collected in the experimental spectra. 

Previous ^tSR work on (Ba,K)Fe2As2 found regions of 
phase-seperated magnetic order, with spontaneous muon 
precession in ZF-/LtSR spectra [17]. In a recent ZF-/^SR 
study of Ba(Feo.gCoo.i)2As2|22], weak, temperature de- 
pendent relaxation was also observed. The presence of 
relaxation due to magnetism largely precludes detailed 
analysis of the properties of the vortex lattice and there- 
fore it is important to check that such relaxation is ab- 
sent. We performed ZF-/xSR measurements at temper- 
atures both above T c and at T= 2 K, finding that the 
spectra were identical, showing only weak temperature 
independent relaxation, such as would be expected for 
nuclear dipole moments. 

Fig. 1 shows a fast Fourier transform (FFT) of the 
TF fiSR signal measured in an applied field of 0.1 T at 
T= 1.6 K. The anisotropic lineshape is characteristic of 
an Abrikosov vortex lattice and indicates the presence 
of at least locally well-ordered vortices within the su- 
perconducting state. The lower cutoff in the lineshape 
corresponds to muons landing at the center of the lat- 
tice (furthest from vortex cores), the maximum in the 
lineshape comes from muons at the saddle point midway 
between vortices while the cutoff at high field comes from 
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FIG. 1: FFT of the TF fxSR signal in B ext =0.1 T at T=1.7 K. 
The anisotropic lineshape is characteristic of a well-ordered 
vortex lattice. 



muons at the vortex cores. The overall width of the linc- 
shape is dominated by the magnetic penetration depth A, 
while the high field cutoff mainly reflects the coherence 
length £. The high field cutoff is most evident in the 
highest applied fields which is furthest from the London 
limit where the field would actually diverge at the vortex 
cores, giving a high field tail with no cutoff. 

The data was analyzed by simultaneously fitting sev- 
eral runs to an analytical Ginzburg-Landau model [TO] for 
an applied field of 0.1 T, to determine the Ginzburg- 
Landau parameter, n = A/£. We found k — 44; this 
value was held fixed for the remainder of the analysis, 
noting that measurements at lower fields are not sensi- 
tive to k as long as it remains large. Small-angle neu- 
tron scattering (SANS) measurements [23J of vortices in 
applied fields above 0.2 T found a highly disordered vor- 
tex arrangement. We included the effects of vortex lat- 
tice disorder in our analysis via an additional gaussian 
broadening of our /iSR spectrum [301 [21], where we as- 
sumed that this broadening was proportional to 1/A 2 as 
observed in previous studies of cuprates and other high 
k superconductors [15] . Our fitted result shows that the 
rms deviation of the vortex positions ((s 2 ) 1 ^ 2 ) relative to 
the vortex separation (Lq) was (s 2 ) 1 / 2 /io w 2.4 % at 
1.7 K in H ex t = 0.1 T. This relatively small amount of 
disorder as seen in the field distribution is somewhat in 
contrast with the SANS results. However, the disorder 
is generally greater in higher fields (as used in the SANS 
study) and may also reflect the fact that //SR, as a real 
space probe is less affected by a loss of true long range 
order than is a reciprocal space probe such as SANS. 
We note that Bitter decoration measurements [23] (in 
low fields) provide evidence for at least regions of well 
ordered vortices. At this point, we were able to fit the 
rest of the parameters in the /iSR signal, including the 
penetration depth, X(H,T). Results of this analysis for 
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1/A 2 oc n s (where n s is the superfluid density) are shown 
in Fig. [2] for applied fields of 0.1 T, 0.05 T and 0.02 T 
where it can be clearly seen that both the temperature 
dependence and zero temperature value n s (T — ► 0) de- 
pend on the applied field. 

The temperature dependence of the superfluid den- 
sity reflects how easily thermal fluctuations are able cre- 
ate quasiparticles. In conventional weak coupling BCS 
theory one finds that the low temperature behaviour 
of n s becomes exponentially flat, whereas the presence 
of gap nodes is reflected in a low temperature power 
law behaviour (such as T-linear in the case of the d- 
wave cuprates). We first fit to a power law, n s (T) = 
n s (0)[l - (T/T c y\. The results of this fit are indicated 
by the dashed lines in Fig. [2] We found that the fitted 
power increases with increasing field, from 1.62±0.03 at 

0. 02 T to 2.27±0.03 at 0.1 T while the superfluid den- 
sity at T = decreased with field, from 0.2047±0.0015 
fim- 3 to 0.1301±0.0005 ^?n~ 3 . We note that this model 
provides a fairly good fit to the data, but the strong field 
dependence of both the superfluid density and the power 
are somewhat surprising for a range of fields much less 
than H c2 jTTj which encourages us to try different possible 
models to characterize the superfluid density. 

ARPES measurements of the iron pnictides have 
demonstrated that the Fermi surface has multiple sheets 
and that there are at least two different superconduct- 
ing gaps on these bands and as a result, a multi-band 
model of the superconducting state has been proposed 
for these systems [1 EZl I2H EHl EH]- Multigap su- 
perconductivity has been seen in a number of systems 
including MgB 2 [31] . which is an s-wave superconductor 
with a Tc of 39K and where theoretical calculations pre- 
dict that the smaller gap is an induced gap associated 
with the 3D bands, and the larger gap is associated with 
the superconducting 2D bands [25]. The possibility of 
other materials exhibiting multi-band superconductivity 
was realized with the discovery of the two-band supercon- 
ductor NbSe 2 by ARPES measurements [53] . This situa- 
tion has been clearly observed in /iSR measurements [3D? 

1. Like MgB2, this material had two gaps of very dif- 
ferent sizes, which reside on two sets of Fermi surface 
sheets [33], Multi-band superconductivity has been pro- 
posed in many other materials, including CeCoIn5[35, 
135] , BaNi 2 P 2 [37J and PrOs 4 Sb 12 [35] . 

We next fit the superfluid density to a phenomenolog- 
ical two-gap model [39, 40J which has been employed in 
a previous fiSR study of LaFeAs(0,F), Ca(Fe,Co)AsO, 
and fBa.K)FeoAso.[4"T] 



n s (T) = n s (0)-w8n s (A 1 ,T)-(l-w)-Sn 3 (A 2 ,T) (1) 

where w is the relative weight for the first gap, Ai. 
Here, the gap functions are given by, 
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FIG. 2: (Color online) Combined two-gap fit (solid line) and 
power law fit (dashed line) to the superfluid density measured 
from the TF jj.SK measurements at 0.02 T (black diamonds), 
0.05 T (red crosses) and 0.1 T (green squares). 



Sn(A, T) = [ X f(e, T) ■ [1 - f(e, T)]de (2) 

kbT J 

where f(e,T) is the Fermi distribution given by, 



/( e ,T) = (l + e ^ 2 +A(T)Vfe B T ) -l 



(3) 



Here, A* (i — 1 and 2) are the energy gaps at T = 0, 
and Aj(T) were taken to follow the standard BCS tem- 
perature dependence. The size of the gaps, Ai and A 2 , 
and T c were fit globally, while n s (0) and the weighting 
factor, w, were allowed to be field-dependent. 

The results of this analysis are shown by the solid lines 
in Fig. [2] It can be seen that this gives a very good fit, 
with a x 2 that are approximately half of that for the 
power law fit. Based on other experimental support for 
a two-gap model, this gives support to our choice of fit . 

From the fit, we obtained the values of the gaps 
2A 1 /k B T c = 3.768 and 2A 2 /k B T c = 1.565, and T c = 
22.1 ± Q.2K . The larger of the two gap values is quite 
close to the BCS value, while the smaller gap is roughly 
half the BCS gap. These gaps are lower than has been re- 
ported for other iron pnictide compounds, which presum- 
ably indicates that the strength of the superconducting 
pairing varies from system to system. 

The relative weighting factor for the larger gap in- 
creases with B app , from w — 0.655(7) at 0.02 T, w = 
0.766(6) at 0.05 T to w = 0.909(4) at 0.1 T. This in- 
dicates that the applied field acts to weaken supercon- 
ductivity on the bands with the smaller gap. The super- 
fluid density at T = 0, n s (Q) , decreases with B app , giving 
n a (Q) = 0.199(1), 0.1441(4), 0.130(1) for these fields re- 
spectively. 
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FIG. 3: (top) Magnetic field dependence of the magnetic 
field penetration depth in Ba(Feo.93Coo.o7)2As2 measured at 
T= 1.7 K. (Bottom) 1/A 2 oc superfluid density, n s , measured 
at T=1.7K. 



In addition to measuring the full temperature de- 
pendence of the penetration depth for 3 applied fields, 
we also studied the magnetic field dependence of n s at 
T=1.7K (following field-cooling through Tc for each ap- 
plied field). We simultaneously fit pairs of /iSR spectra 
at T= 1.7 K and in the normal state using the same 
procedure as before to obtain the penetration depth at 
T= 1.7 K. The results of this analysis are shown in Fig. [3] 
We can see that the superfluid density shows a small peak 
near 0.02 T, then decreases with increasing field, which is 
in agreement with the trend found in our full temperature 
scans. The peak at low fields may reflect the proximity 
of the applied field to the lower critical field, which is 
about 0.007 T at Y= 5 K[T3]. The penetration depth 
tends to a constant value of around 300 nm for higher 
fields. Recalling the field dependence of the relative gap 
weighting factor w, this high field saturation behaviour 
presumably reflects the loss of superconductivity on the 
bands with the smaller gap. 

In conclusion, we have performed transverse field 
/iSR measurements in the mixed state of single crystal 
Ba(Feo.93Coo.o7)2As2. Employing an analytic Ginzburg- 
Landau theory to obtain the penetration depth A and 
Ginzburg-Landau parameter k we find that the super- 
fluid density n s oc 1/A 2 can be well-described by a 
s-wave two gap model where the corresponding field- 
independent gaps are 2Ai = 3.7ksT c and 2A2 = 1.6fcsT c 
respectively. 

We thank Jeff Sonier for sharing his computer pro- 
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